Vol. 13 No.1 天 文 研究 与 技术 Chaa Gel ERAT 


Jan., 2016 ASTRONOMICAL RESEARCH AND TECHNOLOGY 2016 4F 1 H 


CN 53-1189/P ISSN 1672-7673 


AKA MP-CLEAN 方法 对 扩展 源 图 像 重 建 的 比较 


ake 1,2 1,2 
TRIG AR? | GEE 
(1. 清华 大 学 工程 物理 系 天 体 物 理 中 心 ， 北 京 100084; 2. 粒子 技术 与 辐射 成 像 
教育 部 重点 实验 室 (清华 大 学 )， 北 京 ”100084) 


摘要 : 在 射电 综合 成 像 领 域 ， 通 常 需 要 用 退 卷 积 的 方法 补 全 频 域 稀疏 的 采样 。 由 于 扩展 
源 的 频 域 信息 更 为 丰富 ， 要 补 全 这 些 信息 相对 于 点 源 来 说 更 难 ， 因 此 扩展 源 的 图 像 重 建 是 射 
电 综合 成 像 领域 的 一 大 难点 。 为 了 探索 射电 干涉 扩展 源 图 像 重 建 方法 的 特点 ， 将 最 大 类 与 一 
种 加 速 的 CLEAN 方法 (文中 称 之 为 Multi-Point CLEAN, MP-CLEAN) 对 扩展 源 的 干涉 阵 模 拟 
数据 的 图 像 重建 进行 了 比较 。 通 过 比较 ， 发 现 对 于 同样 的 观测 数据 ， 两 种 方法 都 能 较 好 地 重 
建 图 像 ， 但 MP-CLEAN 方法 的 旁 闪 祛除 效果 与 重建 效果 优 于 最 大 类 方 法 ， 而 且 在 模拟 数据 
重建 中 MP-CLEAN 方法 的 总 体 速 度 比 最 大 快 3 倍 以 上 。 最 后 ， 在 讨论 部 分 通过 研究 两 种 
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射电 天 文 是 天 文学 的 一 个 分 文 ， 射 电 天 文 顾 名 思 义 是 通过 探测 射电 波 观测 天 体 ， 从 而 研究 天 体 的 
一 门 科学 。 射 电波 的 波长 相对 于 可 见 光 长 得 多 ， 这 使 得 射电 观测 的 角 分 辩 率 受到 很 大 的 限制 ， 以 
6 cm 的 射电 波 为 例 ， 其 波长 是 可 见 光 (波长 400 至 760 nm) 的 十 万 倍 左 右 ， 也 就 是 说 ,一 个 1 m 直径 
的 光学 望远镜 的 分 辨 率 相当 于 直径 约 100 km 的 射电 望远镜 ， 显然 分 辨 率 是 射电 天 文 发 展 的 一 个 很 大 
的 阻碍 。 从 19 世纪 40 年 代 起 ， 干 涉 成 像 技 术 的 发 展 克 服 了 这 种 局 限 。 干 涉 成 像 技术 的 创始 人 Martin 
Ryle 也 因此 获得 了 1974 年 的 诺 贝 尔 奖 。 接 下 来 的 数 十 年 ， 一 些 大 型 的 射电 干涉 阵列 (如 One-Mile 
Telescope, Green Bank Interferometer, Westerbork Synthesis Radio Telescope, Very Large Array ) 陆续 建成 ， 
大 大 推动 了 射电 天 文学 的 发 展 。 如 今 的 射电 干涉 望远镜 已 经 路 入 世界 上 分 辨 能 力 最 高 的 望远镜 行列 ， 
例如 阿 塔 卡 玛 上 毫米 / 亚 上 毫米 波 阵 列 望 远 镜 ( Atacama Large Millimeter Array, ALMA) 的 空间 分 辨 能 力 已 经 
是 哈 勃 空间 望远镜 (Hubble Space Telescope ) 的 5 倍 左右 。 

干涉 阵 是 利用 干涉 原理 得 到 可 见 度 数据 (Visibility Data) ， 这 里 的 可 见 度 数据 可 等 效 理 解 为 图 像 的 
空间 频谱 数据 。 由 于 干涉 是 来 自 两 个 信号 的 相干 ， 所 以 得 到 的 可 见 度数 据 与 观测 天 线 的 数量 及 相对 位 
置 分 布 有 关 。 受 材料 、 建 造 技术 及 资金 等 限制 ,现实 中 只 能 建造 天 线 数目 和 尺寸 都 有 限 的 干涉 阵 。 但 
是 ， 随 着 干涉 技术 的 不 断 发 展 和 完善 ， 人 们 已 经 开始 筹建 越 来 越 大 的 阵列 ， 例 如 我 国 的 “天 籁 计划 ” 
等 ， 由 此 也 带 来 了 数据 处 理 的 挑战 1。 干涉 阵 的 排 布 对 成 像 有 较 显著 的 影响 ， 这 方面 国内 已 有 所 研 
究 2。 由 于 射电 王 涉 阵 的 观测 数据 是 对 图 像 频 域 空间 的 稀 玻 采样 ， 因 此 如 果 直 接 对 获取 的 可 见 度 数 
据 进 行 依 里 叶 逆 变换 ， 得 到 的 图 像 往往 很 模糊 ， 通 常 存在 很 强 的 旁 淮 。 为 了 重建 清晰 的 图 像 ， 通 常 需 
要 进行 反 卷 积 (Deconvolution) 。 在 射电 天 文 图 像 重 建 领域 ， 常 用 的 几 种 反 卷 积 方法 有 CLEAN 方法 
及 其 衍生 方法 (如 Multi-scale CLEAN! 、Multi-resolution CLEAN 等 ) MRAKA A KAN ETA 
(如 Multi-scale MEM “等 )。 射 电 干 涉 成 像 中 ， 扩 展 源 的 图 像 重 建 是 其 中 的 难点 ， 普 通 的 CLEAN 方法 
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虽然 对 简单 的 点 源 重 建 效果 较 好 ， 但 对 于 复杂 的 大 面积 弥散 源 来 说 重建 时 计算 量 大 ， 收 敛 速度 慢 ， 重 
建 稳定 性 差 ， 得 到 的 结果 通常 不 理想 !'”" 。 具 有 代表 性 的 适用 于 扩展 源 干 涉 阵 数据 的 重建 方法 主要 有 
两 种 : 一 种 是 在 CLEAN 的 基础 上 发 展 起 来 的 多 尺度 方法 ， 以 Cornwell 的 Multi-Scale CLEAN’? 为 代 
表 ， 一 种 是 以 最 大 焙 方 法 (MEM) 为 代表 。 本 文 使 用 一 种 简单 的 加 速 CLEAN 的 方法 ,与 CLEAN 方 
法 每 次 选择 残余 图 中 的 最 亮点 不 同 ， 每 次 选择 大 于 最 亮点 y(0<y<1) 倍 的 多 个 点 ， 而 不 仅仅 是 一 
个  ， 这 种 方法 本 文 称 为 Multi-Point CLEAN 方法 (简称 MP-CLEAN 方法 ) 。 这 种 方法 既 加 速 了 
CLEAN， 又 在 一 定 程 度 上 利用 了 相 邻 像素 值 之 间 的 物理 联系 ， 适 合 进 行 复 杂 扩 展 源 的 重建 。 本 文 之 
所 以 选择 用 最 大 信 方 法 与 MP-CLEAN 方法 进行 对 比 ， 是 由 于 MP-CLEAN 方法 的 简单 性 ， 它 并 不 属于 
多 尺度 方法 ， 同 样 地 ， 最 大 入 方 法 也 不 属于 多 尺度 的 方法 ， 二 者 同属 于 扩展 源 重建 的 初级 方法 。 因 此 
可 以 将 这 两 种 比较 初级 的 方法 进行 比较 ， 通 过 对 比 两 种 扩展 源 重 建 方法 ， 可 以 了 解 两 种 方法 的 特点 ， 
并 为 进一步 完善 和 改进 射电 干涉 扩展 源 图 像 重 建 方法 提供 启发 。 


1 模型 的 建立 


为 了 简化 模型 ， 本 文中 认为 天 图 的 可 见 度数 据 与 天 图 的 空间 频谱 等 效 ， 即 直接 使 用 传 里 叶 频 谱 数 
据 作为 可 见 度数 据 。 任 意 给 定 一 个 亮度 分 布 x, y) ， 甚 可见 度 可 以 用 (1) 式 表示 : 


V(u, v) = fra, y) hey) dedy, (1) 


IF, V(u, v) 为 可 见 度数 据 ; Ia, y) 为 亮度 分 布 图 像 ， 即 可 见 度 是 图 像 的 傅 里 叶 变 换 。 射 电 望 远 镜 
的 观测 等 效 于 在 天 图 的 频 域 进行 采样 。 如 果 准 确 获得 可 见 度 的 全 部 采样 ， 可 以 直接 利用 传 里 叶 逆 变换 
RAS T(x, y) ， 如 (2) 式 : 


I(x, y)= Iv, v) e? dudo . (2) 
然而 实际 数据 往往 是 对 可 见 度数 据 的 部 分 采样 ， 采 样 函 数 如 (3) 式 : 
S(u, v) = $,  @,8(u - u) 6(v - %,), (3) 


式 中 ，o, 为 加 权 因 子 ， 在 本 文中 全 部 取 为 常数 1， 即 认为 采样 完全 准确 有 效 ; 6 函数 满足 6(0) =1, 4 
x 关 0 时 ，6(x) =0。 于 是 得 到 可 见 度数 据 : 


Vu, v)=S(u, v) Vu, v) = 2 Vp v,) 5 (4) 


式 中 ,第 个 采样 点 的 位 置 由 采样 函数 决定 ， 采 样 得 到 的 可 见 度 数据 就 是 六 显然 上 述 过 程 是 不 可 首 
的 ， 信 息 有 丢失 ， 这 就 决定 了 后 面 的 方法 只 能 基于 一 定 的 模型 对 真实 信号 进行 估计 ， 而 无 法 准确 无 误 
地 恢复 复杂 的 真 图 。 

在 重建 时 ,已 知 数据 是 采样 得 到 的 可 见 度数 据 和 采样 函数 (在 图 像 域 表 达 就 是 已 知 脏 图 和 点 扩展 
函数 ) ， 而 需要 重建 真实 图 像 I(%, y) 。 

如 果 考 虑 观测 中 存在 的 噪声 ， 可 见 度 观测 数据 可 以 用 下 式 表示 : 


V(u,) = 之 I(x, exp(2aj u, x;) + nm, (5) 


AP, Vu) RIS ILEWO; x) 表示 真实 亮度 分 布 ; wee UV FER) ER k POR 
位 置 向 量 ; x; 是 图 像 平 面 第 i 个 像素 的 位 置 向 量 ; n, 表示 UV 平面 上 第 大 个 采样 点 的 噪声 。 有 了 (5 ) 式 
描述 的 关系 ， 要求 出 真实 亮度 分 布 ， 需 要 用 反 卷 积 方法 。 
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2 重建 方法 


2.1 RAMIE 
FES ABATE BEA HE SL, OC 11 1991 Ep BSR, (6) sh 
H, = E Lb, = m = bla, (6) 
NF, H, 表示 图 像 的 (相对 ) 焙 ; b 表示 待 重建 图 像 ，m 表示 图 像 模型 ，i 表示 第 i 个 像素 。 
利用 这 种 人 的 定义 ,结合 观测 数据 ， 就 能 求解 真实 亮度 分 布 。 在 观测 数据 的 约束 下 求 满足 粹 最 大 
的 解 ， 可 以 利用 拉 格 朗 日 乘 子 法 构造 目标 函数 : 
J=H, -a° - EQ’)], (7) 
式 中 , XP 为 观测 可 见 度 数据 与 重建 图 像 的 可 见 度 数据 各 点 偏差 大 小 的 归 一 化 平方 和 ; EOC) HX? 的 期 
望 值 , X? 具体 表达 式 为 


5 Dee | 
X= ; (8) 


O, 
RP, o, 是 权重 因子 (本 文中 取 为 1) V, 是 观测 获得 的 可 见 度 数据 ; V EEEE RERE 
数据 ; o? 为 观测 获得 数据 噪声 的 方差 。 
使 目标 函数 取得 最 大 值 的 图 像 b 就 是 待 求 的 解 ， 函 数 的 极 值 问题 有 多 种 求解 方法 ， 这 里 不 妨 转 化 
为 求 目标 函数 梯度 的 0 点 ， 即 可 利用 牛顿 法 求解 ， 迭 代步 长 的 具体 表达 式 由 下 式 得 到 5 ， 


aj 1 
Be Ga aay (9) 
A, n 从 0 开始 取 值 ，5(0) 即 为 已 知 给 定 的 先 验 图 像 mw， 本 文 所 做 模拟 中 的 先 验 图 像 m 是 一 张 所 有 
点 取 值 为 1 的 图 像 ; 为 目标 函数 ;a CR; q 取 点 扩展 函数 的 主峰 峰值 。 

实现 步骤 2 : 

(1) 给 出 初始 重建 图 像 5,， 图 像 模 型 m, XES b, 及 m 所 有 像素 的 值 等 于 1; 

(2) 对 重建 图 像 的 傅 里 叶 域 用 采样 函数 进行 采样 ， 得 到 重 采 样 数据 V, 

(3) 利 用 重 采 样 数据 、 观 测 采 样 数据 及 采样 函数 的 有 关 参 数 ， 代 和 (9) 式 ， 求 出 迭代 步 长 Ab; 

(4) 更 新 重建 图 像 . b=b+Ad; 

(5) 判 断 近 代 终 止 条 件 是 否 达 到 ， 即 归 一 化 卡 方 值 开始 小 于 等 于 1。 如 果 达 到 终止 条 件 ， 则 输出 
图 像 b 即 为 重建 图 像 ， 如 果 未 达到 ， 则 返回 第 (2) 步 ，; 

(6) 将 重建 图 像 卷 积 与 最 大 基线 长 对 应 分 辨 率 的 高 斯 函数 即 得 到 最 终 的 最 大 炉 重 建 结 果 。 
2.2 MP-CLEAN 方法 

在 射电 天 文 图 像 重建 中 ，CLEAN 方法 是 最 常用 的 反 卷 积 方法 之 一 ， 文 [10] 于 1980 年 提出 了 一 种 
可 以 同时 对 多 个 点 进行 CLEAN 来 加 速 CLEAN 的 方法 ,为 了 描述 和 比较 方便 ， 本 文 称 之 为 Multi-Point 
CLEAN( MP-CLEAN)。 在 应 用 MP-CLEAN 时 ,需要 一 条 经 验 性 假设 ， 即 有 采样 较 好 的 数据 对 应 点 扩 
散 函 数 的 主流 比 劳 关 高 很 多 ， 即 使 有 旁 兴 的 又 加 ， 不 同 亮 源 同时 作用 在 同一 位 置 的 概率 不 大 ， 形 成 假 
源 的 峰值 也 不 会 太 大 。 这 就 为 同时 选择 多 个 真实 的 源 提供 了 可 能 ， 因 此 可 以 在 每 一 步 迭代 中 选择 大 于 
最 亮 值 yY(0<y 志 1) 倍 的 多 个 点 ， 而 不 仅仅 是 一 个 点 中 来 加 速 CLEAN。 每 一 次 迭代 y 允许 的 取 值 范围 
可 能 不 一 样 , 但 是 ， 文 中 为 了 简便 ，y 在 同一 次 重建 中 都 取 一 个 较为 保守 的 常数 ( 注 y=1 时 为 
CLEAN 方法 ) ， 参 数 的 选取 在 文中 第 4 部 分 详细 讨论 。 下 面 是 MP-CLEAN 方法 的 实现 步骤 1. 

(1) 给 出 一 个 初始 模型 M,， 并 令 M 所 有 元 素 为 0; 

(2) 利 用 采样 函数 求 出 由 模型 得 到 的 可 见 度数 据 ; 
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(3) 将 模型 可 见 度 数据 与 观测 数据 作 差 ， 再 做 傅 里 叶 逆 变 换 ， 取 实 部 得 到 残 差 图 R; 

(4) 找 出 残 差 图 RR 中 的 最 大 值 WM4X; 

(5) 找 出 残 差 图 中 大 于 MAX y(0<y 和 1) 倍 的 点 ， 作 为 MP 图 (与 模型 尺寸 一 样 的 图 ， 图 中 的 点 除 
了 找 出 的 这 些 点 与 原 图 位 置 及 像素 值 一 样 ， 其 余 像素 点 全 为 0) ; 

(6) 更 新 模型 : 在原 模型 上 加 MP 图 的 B(0<8<1) 售 ; 

(7) 判断 是 否 达 到 人 迭代 终止 条 件 : 残 差 图 对 应 的 可 见 度数 据 的 均 方 根 小 于 等 于 估计 噪声 水 平 ( 即 
卡 方 小 于 等 于 1 时 终止 迭代 ) 。 如 果 未 达到 终止 条 件 ， 则 返回 第 (2) 步 。 

(8) 将 结果 卷 积 与 最 大 基线 长 对 应 分 辨 率 的 高 斯 函数 即 得 到 MP-CLEAN 重建 结 


3 数值 模拟 及 分 析 


3.1 模拟 数据 的 生成 

模拟 数据 的 生成 从 干涉 阵 的 排列 及 组 合 开 始 ， 采 样 函数 由 天 线 阵列 及 观测 等 因素 确定 ， 本 文 模拟 
所 用 的 采样 函数 由 排列 在 同一 条 直线 上 的 天 线 组 成 ， 其 位 置 参考 韦 斯 特 博克 综合 孔径 射电 望远镜 
(Westerbork Synthesis Radio Telescope, WSRT) 的 天线 位 置 设置 。 观 测 时 利用 地 球 自身 的 旋转 ， 一 条 基 
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线 对 应 的 点 在 采样 空间 扫 过 一 条 弧 线 。 为 了 简便 ， 用 一 条 弧 线 表示 这 种 观测 。 观 测 时 间 越 长 ， 扫 过 弧 
有 对 称 性 ， 因 此 扫描 过 程 等 效 于 绕 着 中 心 直 径 上 的 点 旋转 采样 ， 这 里 扫描 的 弧 线 是 椭圆 孤 ， 椭 圆 的 长 
轴 与 短 轴 之 比 为 1:sing (这 里 的 8 为 被 观测 源 所 在 的 赤 纬 ， 文 中 0 取 60° 为 例 ) 。 
尺寸 为 1200 x 1200， 这 里 每 个 像素 大 小 设 为 1 Original image original image 
arcsec X 1 arcsec ) 做 传 里 时 变换 ， 即 得 到 可 见 度 | |° 
对 可 见 度 采 样 并 加 上 噪声 就 可 以 得 到 模拟 的 可 见 
度 观测 数据 (图 1 右 下 图 )。 干涉 阵 的 最 大 基线 
所 选 的 观测 波长 A 约 为 6cm， 得 到 的 分 辨 率 约 RA/arcsec 
为 4.4 arcsec X 5. 1 arcsec。 由 于 观测 数据 分 辨 率 
真 图 ( 真 图 与 原 图 不 同 ， 真 图 是 原 图 卷 积 高 斯 函 
数 得 到 的 分 辨 率 为 4.4 arcsec x 5. 1 arcsec) 的 比 
由 于 模拟 可 见 度数 据 是 复数 ， 因 此 选择 的 噪 
声 也 必须 取 复 数 形式 ， 其 幅度 取 服 从 均 方差 为 5 
对 值 ， 其 相位 取 [0， 27] 上 的 均匀 分 布 。 在 后 面 Al 左上 : 模拟 选用 的 扩展 源 ， 单 位 : Jy; 左下 : 12h 观 
的 模拟 中 假设 知道 噪声 的 大 小 ( 即 预 估 的 噪声 水 测 对 应 的 采样 函数 ， 又 称 为 UV Žž (uv-coverage), 


线 长 度 越 长 ， 考 虑 到 地 球 自转 的 周期 性 ，12 小 时 等 效 180 "的 扫描 。 由 于 实 信号 对 应 的 可 见 度数 据 具 
选 一 张 扩展 源 作为 原 图 (图 1 左上 图 : AR je 
(图 1 右上 图 )， 再 利用 采样 函数 (图 1 左下 图 ) 
长 度 决定 了 其 分 辩 率 (在 波长 一 定 的 情况 下 )， a 
的 限制 ， 为 了 统一 分 辨 率 ， 在 后 面 的 重建 图 像 与 
较 中 ， 统 一 卷 积 一 个 归 一 化 的 二 维 高 斯 函数 。 
(单位 : Jansky, Jy)、 均 值 为 0 的 高 斯 分 布 的 绝 
两 坐标 轴 单 位 都 是 k 和 (1kA 二 60 m); GE: RA 


a l 可 见 度数 据 的 幅度 的 对 数 ， 右 下 : 模拟 观测 ( 采 
由 于 采样 函数 没有 对 可 见 度 的 原点 (直流 分 样 ) 得 到 的 可 见 度数 据 ( 已 加 噪声 ) 的 幅度 的 对 数 
量 ) 进行 采样 ’ 即 总 流量 需要 估计 > 而 这 里 模拟 Fig.1 Top-left panel; extended source chosen for this paper 
中 流量 估计 认为 准确 ， 即 在 频 域 采样 时 对 可 见 度 (unit: Jy); Bottom-left panel; the UV coverage of a 
的 原点 多 采 了 一 点 。 数 值 模拟 中 输入 信息 是 图 1 12h observation; Top-right panel: the logarithm of 
右 下 图 对 应 的 采样 得 到 的 可 见 度数 据 ( 即 模拟 观 the visibility amplitude to the original image; Bottom- 


测 数据 ) 、 佑 计 的 噪声 、 估 计 的 流量 及 采样 函数 。 left panel: the logarithm of the sampled visibility 
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ALARA, MY Wok tine aka MP-CLEAN 
等 方法 对 真 图 进行 重建 。 然 而 由 于 采样 数据 本 身 
的 缺失 ， 对 于 模拟 中 用 到 的 复杂 扩展 源 而 言 ， 任 
何方 法 重建 的 结果 都 不 可 能 恢复 到 与 复杂 的 原始 
图 像 完 全 一 致 ， 因 为 任何 反 卷 积 方法 都 可 看 成 在 
缺失 的 频 域 采样 处 插值 ， 重 建 结果 只 能 是 基于 一 
定 的 模型 对 真 图 的 近似 估计 。 
3.2 模拟 重建 比较 

本 文 所 用 模拟 中 a, B, y 3 个 参数 的 选取 : 
Qa=5、B=0.1、y=0.6( 参 数 的 选取 可 参看 第 4 
部 分 的 讨论 ， 选 择 这 些 参数 的 基本 原则 是 使 得 
收敛 的 卡 方 尽 可 能 地 接近 1 且 保 证 重建 的 稳定 
性 ) 。 

EE A JL SB, ET WAP le 
及 MP-CLEAN 对 模拟 的 可 见 度数 据 进 行 重建 ， 
得 到 的 重建 结果 如 图 2。 图 2 展示 了 最 大 粹 方法 
及 MP-CLEAN 方法 利用 8h 和 12 h 的 模拟 观测 
数据 重建 的 结果 。 可 以 发 现 直 接 利用 采样 数据 作 
傅 里 叶 逆 变换 得 到 的 脏 图 (Dirty map) 与 真 图 
(Ture image) 差别 很 大 ， 主 要 体现 在 旁 闪 很 多 ， 
而 且 有 条 纹 状 结构 。MP-CLEAN 方法 和 最 大 炉 方 
法 的 重建 结果 在 不 同 程度 上 抑制 了 旁 兴 ， 正 中 间 
的 面 源 及 它 的 延展 结构 都 恢复 得 较 好 ， 尤 其 是 
MP-CLEAN 方法 ， 重 建 的 结果 很 清楚 ， 同 时 淹没 
在 旁 淮 中 的 点 源 也 得 到 了 较 好 的 恢复 , 但是， 最 
大 炉 方 法 依旧 有 一 部 分 源 没 能 恢复 ， 而 且 存 在 微 
弱 可 见 的 旁 兴 。 从 残余 图 来 看 (残余 图 是 重建 图 
像 频 谱 重 采样 与 观测 得 到 的 可 见 度 数据 做 差 之 后 
傅 里 叶 逆 变 换 得 到 ) | BI EY) BRAS K PE He 
稀 看 见 真 图 中 间 面 源 的 影子 ， 也 可 以 看 到 许多 黑 
色 点 源 的 残迹 ， 而 MP-CLEAN 方法 得 到 的 残余 
图 中 几乎 看 不 到 这 些 踪 迹 ， 这 说 明 MP-CLEAN 
重建 更 胜 一 筹 。 对 比 真 图 、 脏 图 以 及 重建 结 
总 体 来 说 两 种 方法 的 重建 效果 较 好 ， 重 建 结构 比 
较 清楚 ， 比 脏 图 的 效果 要 好 很 多 。 

为 了 更 直观 地 分 析 它 们 之 间 的 差别 ， 在 RA 
=601 arcsec 处 取 一 剖面 (下 文中 所 有 的 剖面 都 是 
此 位 置 处 的 剖面 ) 如 图 3。 表 1 是 从 剖面 的 残 
差 图 ( 此 处 残 差 图 是 指 真 图 减 去 重建 图 像 得 到 的 
图 ， 注 意 与 前 面 的 残余 图 有 区 别 ) 图 3 中 得 到 的 
几 个 统计 量 ( 由 于 随机 噪声 的 存在 ， 对 于 不 同 次 
的 重建 ， 表 1 至 表 3 中 的 所 有 统计 量具 有 一 定 的 
变化 幅度 ， 但 均值 变化 幅度 小 于 0. 5 JJ ， 标 准 差 
变化 幅度 一 般 不 超过 10% ， 由 于 这 种 变化 很 小 ， 
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图 2 MP-CLEAN ARAZ EM Sh A 12h 模拟 观测 的 
SRR, AL: AA; 右上 : 8h 观测 对 应 的 及 
A; 中 间 和 下 面 各 4 幅 子 图 ， 分别 对 应 8h 和 12h 
观测 数据 的 重建 结果 。 重 建 结果 中 的 4 幅 子 图 分 
NA: BL, RAMERBAR;, AL, RARA 
果 的 残余 图 ; AF, MP-CLEAN 重建 结果 ; GT: 


MP-CLEAN 重建 的 残余 图 


The reconstruction results of the 8h and 12h observa- 


Fig. 2 


tions from the MP-CLEAN and the MEM. Top-left 


panel; the True image; Top-right panel; the dirty 


image made from the 8h observation. The two four- 


subplots at the middle and bottom of the image corr- 


espond to the 8h and 12h observation respectively. 
For each four subplots, Middle-left panel; the MEM 
result; Middle-right panel; the residual of the MEM 
result; Bottom-left panel; the MP-CLEAN result; Bottom- 
right panel; the residual of the MP-CLEAN result 
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并 影响 两 种 方法 结果 的 对 比分 析 ， 因此 本 文具 表 1 重建 结果 剖面 残 差 图 的 均值 及 标准 差 

取 某 一 次 重建 结果 进行 分 析 ) ， 均 值 越 接近 0 且 标 Table 1 The mean and standard deviations of profile 

准 差 越 小 ， 表 明 越 接近 真 图 。 最 大 彤 与 MP-CLEAN differences from the reconstruction 

方法 重建 的 图 像 均值 与 脏 图 相 比 小 得 多 , 说 明 两 

种 方法 的 重建 都 保证 了 总 流量 恢复 得 很 好 。 ee ee ee eee 
从 标准 差 来 看 ，MP-CLEAN 的 结果 无 论 是 8h ae fe ee PR 

还 是 12h 观测 中 ， 其 标准 差 都 要 比 最 大 精 小 。 。 名/ ee ae 

从 剖面 图 可 以 看 出 ， 即 使 在 源 较 弱 的 图 像 边 缘 ， 

MP-CLEAN 的 结果 与 真实 图 像 吻合 得 很 好 ， 而 最 Mean/Jy 16.5 0. 15 -0. 52 
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(在 剖面 RM=601 arcsec HAA MAMA, RAMA MP-CLEAN 的 结果 ) 
Fig.3 The profile difference map of reconstruction results of the 8h (Top subplot) and 12h 


(bottom subplot) observation (true image minus other images at RA=601 arcsec ) 


除了 重建 的 效果 ， 再 看 两 种 方法 的 重建 速度 ， 图 4 2 KIA MP-CLEAN 对 12 h 观测 数据 
进行 重建 的 及 残 差 (绝对 值 ) 最 大 值 随 迭 代步 数 及 时 间 变 化 的 关系 。 从 图 4 中 XX 与 迭代 步 数 之 间 的 
关系 可 以 发 现 ，MP-CLEAN 只 需要 不 到 70 步 , 好 即 达 到 1, MAC a 140 多 步 。 从 X? 随时 间 
的 变化 曲线 可 以 发 现 ，MP-CLEAN 只 需 约 15s, MKTG BE 55s 左右 。 结 合 迭 代步 数 及 时 间 ， 可 
以 算出 MP-CLEAN 平均 每 步 运算 所 需 时 间 为 0.2 s ZA, KIA O0. 4s, HIET IL MP-CLEAN 
的 每 一 次 迭代 的 速度 约 为 最 大 炉 方 法 的 两 倍 ， 而 MP-CLEAN RARESA ERKAT EF, 
最 终 导 致 最 大 焙 方 法 所 需 的 重建 时 间 是 MP-CLEAN 的 3 倍 多 (接近 4 倍 ) 。 
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x2 vs iteration x2 vs Time 
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图 4 他 及 残 差 绝 对 值 最 大 值 随和 迭代 步 数 及 时 间 变 化 的 关系 曲线 


N Fig.4 The X? and the maximum absolute residual versus the iteration or time 


为 了 进一步 检验 这 两 种 方法 的 重建 特点 ， 选 择 更 多 的 原 图 进行 模拟 观测 ， 然 后 重建 ， 采 样 时 间 为 
12 h， 采样 阵列 保持 不 变 。 为 了 体现 与 上 面 的 扩展 源 不 一 样 ， 选 用 的 源 一 个 是 规则 的 源 ( 月 球 ) ， 另 一 
E 个 是 两 个 不 规则 星系 的 原 图 。 重 建 得 到 的 结果 如 图 5， 从 图 5 的 重建 结果 不 难 发 现 ， 与 前 面 分 析 的 重 
r= 建 特点 一 致 : 流量 恢复 均 较 好 ， 都 能 较 好 地 重建 ， 而 且 MP-CLEAN 对 旁 兴 的 祛除 更 好 ， 其 剖面 残 差 
© 的 方差 更 小 。 限 于 篇 幅 原因 ， 关 于 这 些 源 重建 速度 的 曲线 并 未 给 出 ， 这 两 个 源 的 实验 结果 是 : 第 1 个 
W (Source 1) 的 重建 中 ，MP-CLEAN 的 迭代 时 间 不 到 15 s， 而 最 大 炉 则 花 了 近 50 s( 是 MP-CLEAN 的 3 
L); 第 2 个 源 (Source 2) 的 重建 中 ，MP-CLEAN 的 欠 代 时 间 不 到 30 s， 而 最 大 信 的 迭代 时 间 却 超过 
T 150 s(# MP-CLEAN 的 5 倍 多 ) 。 模 拟 结果 表明 MP-CLEAN [hie Ke 3 倍 以 上 ， 与 前 面 的 结果 保 


持 一 致 。 
4 讨 论 


在 3. 2 PEZA TKR MP-CLEAN 的 重建 结果 ， 并 且 进 行 了 详细 的 对 比分 析 。 然 而 ， 最 
KG MP-CLEAN 方法 的 重建 效果 与 重建 参数 的 选择 有 关 ， 下 面 是 利用 前 面 8 h 的 观测 数据 在 不 同 参 
数 选择 情况 下 的 重建 ， 通 过 对 比 和 分 析 重 建 结 果 ， 进 一 步 了 解 两 种 重建 方法 的 特点 。 

AFC ARK TT PHB oo, 48K 6 中 不 同 a 情况 下 % 的 变化 曲线 ， 可 以 看 到 如 果 不 断 地 
Ka, XX 会 迅速 地 降 到 1， 重建 结果 的 标准 差 较 大 ， 而 如 果 不 断 减 小 a 到 0 附近， 那么 重建 结果 与 
真 图 就 有 较 大 的 整体 偏差 。 这 可 以 从 表 2 中 看 出 明显 的 规律 : 随 着 a 的 增 大 ， 均 值 的 变化 趋向 0， 然 
后 稳定 在 0 附近 ; 从 标准 差 来 看 ， 随 着 a 的 增 大 ， 标 准 差 逐渐 增 大 。 
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图 5 RAMA MP-CLEAN 方法 对 于 不 同 扩展 源 的 重建 结果 。 中 心 线 左 边 为 第 1 个 源 的 重建 结果 及 在 RA=377 arcsec 
处 的 剖面 残 差 曲线 ; 中 心 线 右边 为 第 2 个 源 的 重建 结果 及 在 Dec=601 arcsec 处 的 剖面 残 差 曲 线 


Fig.5 The reconstruction results from the MEM and MP-CLEAN method for different sources. The left panel is the reconstruction 
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results of Source 1 and the plot of the profile difference map at RA=377 arcsec; the right panel is the reconstruction 


results of Source 2 and the plot of the profile difference map at Dec=601 arcsec 


这 与 a 作用 在 观测 数据 上 有 关 ， 随 着 a 增 大 ， 观 测 数 据 这 个 约束 的 权重 越 大 ， 所 以 均值 会 收敛 


到 0， 而 焙 的 权重 就 相对 变 小 ， 业 本 身 是 一 种 平滑 性 约束 ， 业 权 重 减 小 ， 意 味 着 重建 结果 的 振荡 剧 
烈 ， 所 以 标准 差 逐渐 变 大 。 在 a=0.2 i, PERMASI RMT, BELEAK 6 PX? 的 收敛 曲 
线 ， 还 是 从 表 2 中 的 均值 来 看 ， 重 建 结果 与 真 图 有 整体 上 的 偏离 。 从 图 6 中 X 的 变化 曲线 可 以 看 到 ， 
a=0. 2 时 xX 收敛 的 值 远大 于 1， 从 表 2 中 的 标准 差 来 看 ，a=0.2 时 的 结果 振荡 最 小 ， 恰 恰 是 箭 作 为 
一 种 平滑 性 约束 的 体现 。 重 建 结 果 的 剖面 残 差 图 如 图 7， 在 图 7 中 ,除了 a=0.2 时 重建 结果 与 真 图 有 
整体 的 偏 移 ( 从 表 2 中 对 应 的 均值 可 以 看 出 ) ， 难 以 看 出 重建 结果 之 间 较 明显 的 差别 ， 这 也 说 明 最 大 
Wi] HE EZ AR YS NH EA BORE AS PEE 
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下 面 再 看 MP-CLEAN 方法 在 y(0<y<1) 取 不 同 值 时 的 重建 情况 ， 从 图 8 中 的 和 2 随和 迭代 的 变化 曲 
线 可 以 看 到 ，y 越 大 则 MP-CLEAN 方法 对 扩展 源 的 重建 速度 越 慢 ， 图 像 重建 结果 中 的 振荡 也 变 明显 ; 
y 越 小 则 重建 结果 越 平滑 ， 但 是 重建 图 像 中 存在 假 源 的 可 能 性 越 大 ，y 越 小 假 源 的 数目 和 强度 的 期 望 
值 都 会 变 多 和 变 大 ，y 太 小 (如 图 8 中 y=0.1) XX 难以 降下 来 ， 只 能 收敛 到 大 于 1 的 值 。 

重建 结果 的 剖面 残 差 图 如 图 9， 可 以 发 现 y 取 不 同 的 值 对 重建 结果 有 较 显 著 的 影响 ， 对 比 图 9 与 
图 7 RRM, BORMOTIETE a 取 不 同 值 时 ， 重 建 结果 之 间 的 差别 并 不 太 大 ， 而 从 表 2 与 表 3 中 的 标 
准 差 对 比 中 可 以 定量 地 说 明 这 一 点 ， 因 此 MP-CLEAN 方法 对 影响 其 重建 结果 的 参数 y 的 敏感 程度 较 
H, MEREKAM EREE., ME 3 中 可 以 发 现 MP-CLEAN 在 不 同 参 数 取 值 下 均值 都 接近 0, 但 
是 方差 则 呈现 中 间 (y 取 0.5~0.7) 低 ， 两 侧 高 的 总 体 趋势 ， 这 说 明 MP-CLEAN 的 参数 选择 受 两 方面 
的 约束 ， 一 个 约束 是 重建 结果 中 的 假 源 应 当 数 目 少 而 且 强 度 弱 ， 另 一 个 约束 是 结果 的 稳定 性 。 前 一 个 
约束 要 求 y 尽 可 能 地 接近 1， 因 为 这 样 重建 结果 含有 假 源 的 可 能 性 才 会 变 小 ， 后 一 个 约束 则 要 求 y 尽 
可 能 接近 0， 在 这 两 个 约束 的 相互 作用 下 ， 最 后 就 呈现 了 表 3 中 的 方差 中 间 低 、 两 侧 高 的 现象 。 采 样 
越 多 越 全 ，y 可 以 选取 的 范围 越 大 ， 越 有 利于 MP-CLEAN 的 重建 。 如 果 采 样 太 少 ，y 选择 的 区 间 就 会 
越 小 ， 甚 至 可 能 找 不 到 合适 的 y， 此 时 再 使 用 MP-CLEAN 方法 进行 复杂 扩展 源 的 重建 可 能 难以 得 到 好 

的 结果 。 


In (x2) vs iteration In (2) vs iteration 
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Iteration Iteration 
图 6 AB a BAX? 的 自然 对 数 图 8 RE y 取 值 下 X? 的 自然 对 数 
随和 迭代 变化 的 关系 曲线 随和 迭代 变化 的 关系 曲线 


Fig.8 The natural logarithm ofA? versus the 


Fig.6 The natural logarithm of X° versus the 


iteration when @ varies iteration when y varies 
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重建 结果 剖面 的 残 差 图 重建 结果 剖面 的 残 差 图 


Fig.7 The profile differences of the MEM reconstruction Fig.9 The profile differences of the MP-CLEAN recon- 


results within the 8h when a varies struction results within the 8h when y varies 
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R2 RAAH RENTE e 取 值 下 重建 图 像 的 表 3 MP-CLEAN 在 不 同 y 取 值 下 重建 图 像 
剖面 残 差 的 均值 及 标准 差 的 剖面 残 差 的 均值 及 标准 差 
Table 2 The mean and standard deviation of profile Table 3 The mean and standard deviation of profile diff- 
differences using MEM when a varies erences using the MP-CLEAN when y varies 
a Mean/Jy a/Jy y Mean/Jy o/Jy 
0.2 11.1 16.3 0.1 -0.51 20. 0 
0.8 2.37 17.3 0.3 -1.25 21.3 
3.2 0.17 18. 2 0.5 0. 19 16.4 
12.8 -0. 15 18.8 0.7 -0. 20 17.4 
51.2 —0. 26 18.8 0.9 一 0. 97 43.2 


由 于 B 的 取 值 对 MP-CLEAN 的 影响 主要 体现 在 速度 上 ， 对 重建 图 像 效 果 的 影响 不 大 ( 当然 6 很 大 
时 也 可 能 出 现 一 定 的 不 稳定 性 ) ， 因 此 6 对 MP-CLEAN 重建 的 影响 在 此 不 再 歼 述 。 


5 总 结 


/JUN 


AS SCA FARA MP-CLEAN 方法 对 模拟 扩展 源 的 可 见 度 数据 进行 重建 ， 可 以 发 现 : 

(1) 最 大 炉 方 法 有 着 比 MP-CLEAN 更 严格 的 数学 迭代 形式 ，MP-CLEAN 方法 则 带 有 一 定 的 经 验 性 ; 

(2) 对 于 同样 的 采样 数据 ，MP-CLEAN 对 扩展 源 的 重建 比 最 大 箭 方法 重建 效果 好 (使 得 邓 收敛 到 
1 附近 时 各 参数 的 取 值 ) ， 主 要 体现 在 旁 淤 降 得 更 低 ， 重 建 图 像 与 真 图 之 间 差 别 更 小 ; 

(3)MP-CLEAN 方法 的 重建 速度 比 最 大 粹 快 ， 从 模拟 结果 来 看 ， 总 体 上 所 花 时 间 不 到 最 大 粹 的 三 
分 之 一 ; 

(4) 最 大 彤 方法 及 MP-CLEAN 方法 的 重建 效果 与 参数 的 选择 有 关 ， 最 大 粹 方法 对 参数 的 依赖 性 要 
小 于 MP-CLEAN 方法 ， 即 最 大 丧 方法 的 鲁 棒 性 优 于 MP-CLEAN 方法 。 

总 而 言 之 ， 就 本 文 所 用 的 模拟 数据 及 重建 结果 而 言 ， 综 合 模拟 重建 图 像 的 效果 及 重建 的 速度 ， 
MP-CLEAN 比 最 大 粹 方法 有 更 佳 的 重建 效果 及 更 快 的 重建 速度 。 然 而 ， 由 于 MP-CLEAN 对 可 见 度数 
据 采 样 的 要 求 较 高 ， 而 且 需 要 满足 源 的 旁 兴 作用 在 同一 位 置 的 概率 较 小 这 种 经 验 性 的 假设 ， 在 采样 很 
差 ( 点 扩散 函数 旁 欠 峰值 接近 主峰 峰值 时 ) 或 者 在 经 验 性 假设 不 成 立时 (例如 亮 源 的 某 些 对 称 性 分 布 导 
致 旁 闪 全 加 很 严重 ) ，MP-CLEAN 方法 从 原理 上 已 经 不 能 适用 ， 而 最 大 焙 方 法 理论 上 依旧 能 够 用 于 重 
E, DERKI AEAEE MP-CLEAN 好 的 体现 之 一 。 但 MP-CLEAN 方法 在 一 定 程 度 上 挖掘 
T CLEAN 方法 在 同一 尺度 上 的 潜能 ， 可 以 为 完善 和 改进 可 见 度 采样 较 好 的 扩展 源 图 像 重 建 方 法 提供 
更 好 、 更 快捷 的 途径 。 
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Abstract: In the radio synthesis imaging area, the deconvolution is used to fill up the un-sampled points 


at the Fourier domain. The algorithms like CLEAN and MEM often play an important role in the deconvolution. 


As extended sources have more complex visibility, the reconstruction for extended sourcesis difficult. In this 


paper we make a comparison of the MEM ( Maximum Entropy Method) with the speed up CLEAN method 


( named Multi-Point CLEAN in this paper, or MP-CLEAN for short) for extended source reconstruction based 


on the simulated visibility data. The simulated data are produced by sampling the Flourier Frequency domain 


for a true image with the sample function generated from a simulated interferometer array. By analyzing and 


comparing, it is easy to find that both of these two methods work well on the image reconstruction while the 
MP-CLEAN is better than the MEM for getting rid of the side lobes. The MP-CLEAN is at least 3 times faster 


than the MEM in our simulations. A deeper analyzing of the selection of parameters shows that the MEM is 


more independent on the parameters than the MP-CLEAN, which indicates that the MEM is more robust. 


Key words: Image reconstruction; Radio interferometer array; Multi-Point CLEAN; Maximum Entropy Method; 


Extended source 


